Maximum information gain in weak or continuous measurements of qdits: 

complementarity is not enough 
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To maximize average information gain for a classical measurement, all outcomes of an observation 
must be equally likely. The condition of equally likely outcomes may be enforced in quantum 
theory by ensuring that one's state p is maximally different, or complementary, to the measured 
observable. This requires the ability to perform unitary operations on the state, conditioned on the 
results of prior measurements. We consider the case of measurement of a component of angular 
momentum for a qudit (a D-dimensional system, with D = 2 J + 1 where J is the total angular 
momentum). For weak or continuous-in-time (i.e. repeated weak) measurements, we show that 
the complementarity condition ensures an average improvement, in the rate of purification, of 
only 2. However, we show that by choosing the locally optimal control protocol of this type, one 
can attain the best possible scaling, 0{D^), for the average improvement. For this protocol the 
acquisition of information is nearly deterministic. Finally we contrast these results with those for 
complementarity-based protocols in a register of qbits. 
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I. INTRODUCTION 

In the classical world predictability is associated with 
intimate knowledge of a system while unpredictability 
implies surprise. In information theory, the surprisal of 
an outcome k is defined as I{k) = — log(Pfc), were is 
the probability of that outcome The suprisal is an 
important quantity as it quantifies the amount of infor- 
mation one learns from an outcome of a measurement. 
When the average surprisal of an experiment is maxi- 
mized the observer maximizes the extraction of informa- 
tion from the experiment. 

In quantum theory the relationship between pre- 
dictability and knowledge is not clear cut. Even when an 
observer has maximal knowledge about a quantum sys- 
tem (that is, the observer's state of knowledge is pure), 
measurement outcomes may be unpredictable. This is 
the case if measurements on the system are performed in 
a basis which is not the eigcnbasis of the state. (The 
eigenbasis means the basis in which the operator un- 
der consideration, e.g. p, only has entries on the diago- 
nal.) Nevertheless pure quantum states are the most pre- 
dictable states available. It is for this reason, and others 
eloquently explained in Ref. 0, that we take the impu- 
rity to be a measure of information. These considerations 
have lead to information [H and control Q theoretic for- 
mulations of information-disturbance relations. 

Inspired by this work, Jacobs asked the question; how 
quickly on average can information be extracted from 
a quantum system using continuous quantum measure- 
ments in an complementary basis (Other ways of 
quantifyingthe rate of information gain can also be con- 
sidered dilg].) He found, for a widely applicable contin- 
uous measurement model, that measuring in an comple- 



mentary basis led to a speed-up in the purification by a 
factor of two for a two level system Q • (Practically, such 
a rapid purification protocol can be a primitive for a rapid 
state preparation protocol 0, 0]-) Later this was gen- 
eralized to complementary bases for ZJ-dimensional sys- 
tems. In Ref. Q the lower bound on the speed-up for D- 
dimensional systems was found to be Slb = (2/3)(£' + 1) 
when the monitoring was performed in an complemen- 
tary basis. Recently in Ref. [1| the upper bound S\jb 
was proven to be bounded above by for measure- 

ment in any complementary basis. 



In this paper we give an explicit construction of an 
complementary measurement protocol for a qudit (a D- 
dimcnsional system, with D = 2J + 1 where J is the 
total angular momentum) based on the quantum Fourier 
transform (QFT) which achieves a speed-up of Q.2D^. 
Further we show that transforming the state's eigenbasis 
so that it is complementary with respect to the measure- 
ment eigenbasis is not sufficient to achieve this speed-up; 
the speed-up can be as low as 2. In addition to this we 
show, numerically, that unbiased bases can achive the 
upper bound, saturating the inequality S < -D^/2, for D 
even. To define these speed-ups rigorously, it is neces- 
sary to derive rigorous bounds on information extraction 
rate for commuting measurements without control, so we 
also do that in this paper. (For the rate of commuting 
measurements with control see Refs. 1,0). Finally we 
consider measurements in commuting and complemen- 
tary bases for a register of qbits. The complementary 
measurement scheme based on the QFT appears only to 
give a speed-up of 2 in the case of a register. 
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II. FURTHER BACKGROUND AND 
STRUCTURE OF THIS PAPER 

In this section we explain the structure of our paper 
by reviewing the relevant prior work on the subject and 
relating it to the work presented in this article. 

In the quantum information context, information- 
disturbance relations were inspired by the epistemic in- 
terpretation of quantum states as applied to quan- 
tum cryptography 0. The idea was that the informa- 
tion gathering and disturbance should be grounded with 
respect to an observer's state of knowledge, not some 
pre-existing property of the system. To formalize this in- 
tuition Fuchs and Jacobs in Ref. Q considered something 
like the following game. Two observer's Alice and Bob 
initially agree that the state of some quantum system is 
Pq, an impure state. At some later time Charlie will per- 
form a measurement along a randomly chosen axis and 
ask Alice and Bob to predict the outcome of the measure- 
ment. After repeating this many times Charlie rewards 
the individual who makes the most correct guesses. Al- 
ice, being clever and resourceful, has the ability to per- 
form one unitary operation followed by a two-outcome 
measurement. Her aim is to refine her state of knowl- 
edge about the quantum system, thus increasing its pre- 
dictability and accordingly the likelihood of receiving the 
money. However, Alice's morals prevent her from inten- 
tionally sabotaging Bob for the pecuniary reward. Thus 
Alice would like to increase the accuracy of her predic- 
tions without affecting those of Bob. 

This formulation suggests (intuitively) that weak mea- 
surements are the right way to probe the system. From 
the point of view of the weak measurement formalism, 
a projective measurement corresponds to a measurement 
of infinite strength . For finite strength measurements 
on a single qbit it was found that there is a non-trivial 
tradeoff between Alice's information gathering actions 
(the refinement of her predictability) and the disturbance 
to Bob's predictability (his state of knowledge) . 

For a fixed measurement rate the interval of the non- 
trivial tradeoff may be characterized by the angle be- 
tween the Bloch vector and the measurement axis. There 
are two extremes of this tradeoff. Minimal disturbance 
and minimal average information gain occurs when Al- 
ice's axis is aligned with the Bloch vector, that is the mea- 
surement eigenbasis and the state's eigenbasis commute. 
Maximal disturbance and maximal average information 
gain, is attained when Alice's axis is orthogonal to the 
Bloch vector. This strategy corresponds to making the 
outcomes equi-likcly, that is maximizing the suprisal [ij 
of the measurement result. In this situation the bases are 
said to be complementary, maximally non-commuting, or 
unbiased. Thus we conclude that quantum theory does 
not allow Alice to succeed. The ability to measure in dif- 
ferent bases is an essential feature of quantum mechanics 
which provides the richness to the problem considered 
here. 

It is easy to understand how non-commuting measure- 



ments provide maximal information gain and maximal 
disturbance through the following picture. Consider a 
mixed qbit state with its Bloch vector aligned along the 
+ \z) direction. By measuring in the \x) basis the Bloch 
vector will elongate regardless of the measurement out- 
come. This point is illustrated in Fig. [1] (a). Furthermore 
both outcomes are equally likely, which means the aver- 
age information extraction is maximal. Unfortunately 
the post-measurement state is not pointing in the same 
direction as the pre-measurement state, so it will not 
commute with the original state. Conversely if the mea- 
surement basis and the state are along the \z) direction 
then, due to the stochastic nature of the measurement 
outcomes, the average elongation of the Bloch vector will 
be minimal 0, [Tl| . In any particular measurement the 
impurity may increase or decrease. But the disturbance 
to the state will be minimal Q . In Fig. [T] (b) we illustrate 
that one of the possible measurement outcomes causes 
the Bloch vector to shrink whist the other elongates it. 




FIG. 1: A Bloch sphere depiction of the extremes of 
information-disturbance tradeoffs. In both figures the dark 
arrow represents the Bloch vector of the pre-measurement 
state. The arrows labeled with M± correspond to the post 
measurement states conditioned on the random outcome of 
a weak measurement. Figure (a) is measurement of a state 
in a complementary or unbiased basis, in this case the mea- 
surement axis is along the x direction. When this occurs the 
Bloch vector always elongates, as indicated by the light grey 
arrow (the plus outcome rotated for easy visual comparison). 
Figure (b) depicts the commuting measurement. In this case 
if the measurement result M_ is obtained the Bloch vector 
decreases in length. 

While information-disturbance relations are interest- 
ing for single measurements the subject has been thor- 
oughly explored. Less exploration has been undertaken 
in continuous measurement setting due to the increased 
difficultly of the analysis. In order to make some progress 
the focus on disturbance has mostly been ignored (for 
good reasons as wc shall see). That is. the focus has 
been on information gathering. This restriction has led 
to an interesting development in the theory of continuous 
quantum measurements 

To demonstrate an improvement in the extraction of 
information the extraction rate of a standard continuous 
measurement must be defined. We follow Jacobs and de- 
fine the standard measurement as a measurement in a 
commuting basis This is sensible for a number of 
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reasons. Firstly measurement will rapidly destroy the 
coherences of any state that does not commute with the 
measured observable [lH . This can be understood as the 
measurement projecting the state to the eigenbasis of the 
observable. To counter this effect an adaptive measure- 
ment or closed loop feedback control is required. The sec- 
ond reason is a continuous measurement in a commuting 
eigenbasis is essentially a classical continuous measure- 
ment. 

The analysis of disturbance and information gain is 
significantly harder for continuous measurements of D- 
dimensional systems [1, [l^. Minimal disturbance still 
occurs when measuring in the same eigenbasis as the 
state (a commuting basis). A full analysis of this sit- 
uation would require characterization of the distribution 
of information extraction rates, allowing calculation of 
the average information extraction, median, worst case 
and best case. In Sec. IIIII we rigorously derive the av- 
erage purification rate for a commuting measurement of 
a ZJ-dimensional system and provide an asymptotic ex- 
pression. Additionally we give exact results for bounds 
on the spread of the distribution of purities, in Sec. [Ill Dl 
which were heuristically derived in previous work [16[. 

Now that we have defined what a standard measure- 
ment is we can explore what it means to increase the rate 
of extraction of information. In the context of control 
theory one might be interested in revealing or acquir- 
ing the value of some parameter; that is the domain of 
adaptive measurements [13, [l3 • But more typical con- 
trol objectives are to drive the evolution of the state or 
to stabilize the state against noise [iH ITsl - Eol . Thus 
disturbing the prior state is not a problem. This why 
disturbance is not usually considered in feedback control. 
What becomes important is the noise introduced by the 
measurement to the system. Doherty, Jacobs and, Jugn- 
man in Ref. showed, for a continuously monitored 
qbit, that the measurement noise is maximal when the 
eigenbases of the state and measurement axis are unbi- 
ased. 

The study of Wiseman, Mancini and Wang in Ref. [1^ , 
explored state stabilisation in a qbit. The relevant result 
from Ref. [l^l to the present discussion relates to the 
difficulty of stabilising a state which is complementary 
to the measurement basis. They showed that this can 
be achieved in the limit of infinite strength "Bayesian" 
feedback, but not in simpler "Markovian" feedback pro- 
cedures. This is a consequence of the increased measure- 
ment noise as predicted by Doherty et al. Q . Later Ja- 
cobs showed how continuously monitoring the state of a 
qbit in a complementary basis (assuming infinite strength 
feedback) leads to a doubling rate of information extrac- 
tion over monitoring in a commuting basis asymptotically 
i- 

Recently a number of studies have examined the use of 
weak (and continuous) measurement in a complementar- 
ity basis in systems of dimension D greater than two. The 
first study showed that the lower bound on the improve- 
ment in the rate of information extraction when measur- 



ing in any complementary is (2/3)(-D-|-l) Q. The second 
study showed this speed-up is upper bounded by a factor 
We summarise these results in Sec. IIV Al In 
Sees. lIVBI to llVDI of this paper we give an explicit con- 
struction of an unbiased basis protocol for an arbitrary 
finite dimensional system based on the quantum Fourier 
transform. We derive upper and lower bounds on the 
speed-up of the QFT protocol. Surprisingly the numer- 
ics show that the protocol saturates the upper bound 
- for even D. For completeness we mention the recent 
work of Ruskov et al. |2l| , where it was shown that mon- 
itoring a qbit in three complementary bases where they 
claimed a speed-up S' ~ 3. We are currently exploring 
this relationship so we will not discuss it further in this 
article. 

Motivated by some of the above investigations Jacobs 
and Lund argued that the task of state stabilisation is 
best accomplished by measurement in an complementary 
basis; in the regime of strong feedback [l^. This is one 
reason to continue investigations into maximal informa- 
tion extraction, foundational issues aside. In Ref. [l^ 
the authors raised two interesting questions: first, are all 
complementary bases are equally good for information 
extraction when D > 4? Second, is measurement in a 
complementary basis is best for information extraction? 

We reexamine these questions in Sec. |V] and Sec. I VII 
We show that complementary, i.e. equi-likelihood con- 
dition, is not sufficient to guarantee maximum informa- 
tion extraction. Our analysis provides an insight into the 
mechanism for information extraction, which is quantum 
mechanical backaction. From this we explain how maxi- 
mum information extraction is achieved by enforcing the 
complementary condition and arranging the eigenvalues 
in phase space so measurement backaction can have max- 
imum effect. We provide further evidence for the claim 
that not all bases are equally good. Because Jacobs and 
Lund did not examine permutations of the unbiased ba- 
sis it was still unclear from their analysis whether all 
unbiased bases are equally powerful for information ex- 
traction. So we also consider the role of permutations in 
Sec. IVII Our results show that even after maximizing in- 
formation extraction over permutations of the unbiased 
basis, not all unbiased bases are equally good at infor- 
mation extraction. 

In Sec. I VII] we examine complementary measurements 
in a register of qbits using the QFT. We start our anal- 
ysis, in Sec. IVII Al with the case of commuting measure- 
ments 0. Then, in Sec. IVII Bl we review the results of 
Ref. where it was shown that the speed-up bounded by 
2ri,(2"- 1)^1 < 5 < 2ri. In Sec. IVII CI we construct a pro- 
tocol based on the QFT and present results which suggest 
the speed-up is at most 2. We conclude in Sec. I Villi 



III. COMMUTING MEASUREMENTS 

Later in this article we will investigate strategies that 
increase the rate at which information can be extracted 
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from quantum systems. Before these strategies are exam- 
ined, we must fully characterize this rate in the absence 
of feedback. This will provide a benchmark to assess the 
performance of the strategies that aim to increase the 
rate of information extraction. 

In this section we restrict our attention to continuous 
measurement (CM) of a single _D-dimensional quantum 
system, also called a qdit, monitored by a single output 
channel. This output channel is monitored in such a way 
that it results in a diffusive quantum trajectory [l^ . Il4| . 
This can be thought of as a sequence of short-duration, 
weak measurements. Following Fuchs and Jacobs ^ our 
notion of information is characterised by the impurity: 
L[p{t)] = 1 — Tr[p(i)]. Accordingly we derive an exact 
expression for the average impurity {L{t)) and a simpler 
approximate expression in Sec. IIII Al In Sec. IIII Bl and 
Sec. IIII CI we examine the validity of our approximate 
expression. To completely characterise the rate of in- 
formation extraction from commuting measurements we 
look at bounding the distribution of purities in Sec. IIII D] 
Then finally we summarise the key messages in Sec. IIII El 



A. Qdit - impurity analysis 

The starting point for the current analysis is similar to 
that presented in Refs. [1, H, H^. We use a widely appli- 
cable model that describes the evolution of the state p of 
a quantum system due to repeated weak measurements 
of a dimesionless observable X. This model is described 
by the stochastic master equation (SME) [l3, \2§t^ 



dp = dtCc [X] p = 2-fdtV[X]p+ ^27 dW n[X]p, (1) 

where V[A]p = ApA^ ~ UA^Ap + pA^A), H[A]p = 
Ap+pA"! -Tr [{A^ + A)p] pJll, and dW is the increment 
of a Wiener noise process |2a^. It should be noted that 
we have moved to a frame that has enabled us to factor 
out the Hamiltonian evolution. The measurement rate, 
7, determines the rate at which information is extracted, 
and thus the rate at which the system is projected onto a 
single eigenstate of X [l^, [27l - [29l ] . (This means for times 
T 3> 7""*^ we may say that we have performed a projective 
measurement of the observable X.) The measurement 
result in a small time interval [t,t + dt) is denoted 



(2) 



We will denote the continuous measurement record ob- 
tained by the observer integrated up until time t as 



dR{t'). 



(3) 



The unnormalised version of Eq. ([TJ is know as the linear 
SME. We denote this linear SME as [H, H^-IH 



dp = dtCc [X] p = 2jdtV[X]p+ dR U [X] p. (4) 

Here ii[A] p = Ap + pA^ , and the bar over p denotes the 
lack of normalization at all but the initial time. We take 
the initial state to be p = p(0) = I/D where I is the 
D X D identity matrix. Because the initial state matrix 
commutes with X at all times this makes obtaining the 
solution to Eq. (U) simple. This solution is 



p{R,t) = eiip{-A-fXH)exp{2,/2-fXR{t))I/D. (5) 



Using the linear trajectory (see Refs. [30|-|32|) solution, 
i. e. Eq. ([5]) , we now explicitly calculate the form of the D- 
dimensional state matrix. We take the observable X to 
be the z component of angular momentum (represented 
by the operator J^), thus D ~ 2J + 1. In matrix form 
the Jz operator is — diag( J, J — 1, . . . , — J + 1, — J)- 
Using this representation of the unnormalized solution 
of the state evolution can be written as matrix elements 
exponentiated: 



g-47(J-l)"tg2y2^(J-l)fl 



... ^e-47(-J)^*e-2v/2^Jfly (6) 



The normalization factor can be written with a non- 
symmetric sum or a symmetric sum 

, 2J 

JV = Tt[p] = 1 ^ e-47(J-™)^te2y2^(j-™)i? 



m=0 

1 

- y 



D 



(8) 



s=-J 



dR{t) = ^A-i{X{t))dt + dW{t), 



For discussions later we also calculate the probability dis- 
tribution of the result R and choose the symmetric sum 
for the normalization 



where dW is the same Wiener noise process that appears 
in Eq.dl]), and {X{t)) = Tr[Xp(t)]. A quantum trajec- 
tory is a continuous-in-time description of the state con- 
ditioned on the measurement result. Thus, to specify a 
trajectory over the interval [0,r) we must be given the 
measurement results over that interval. Then the Ito def- 
inition pit + dt) — p{t) + dp(t) may be used recursively to 
generate a trajectory. Clearly a single weak measurement 
is described by one application of this relation. 



dRV{R) = dRTT[p{R,t)]P{R) 

J 



dR—^ y , 



-{B.~2^/2^sty /2t 



(9) 
(10) 



Here P{R) has been called the ostensible distribution for 
i? [ilEi. It is given by P{R,t) = cxp {-R"^ /2t)/V^■ 
B^y choosing the symmetric summation for the normal- 
ization one can see that the probability distribution for R 
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contains D peaks centred around the values oc s, where 
s G [J, J— 1, . . . , — J + 1, — J]. For times t » 7~\ the dis- 
tribution is sharply peaked about th ese D values. Each 
of these peaks has a FWHM of 2^2 ln2v^ « 2.35^/1. 
This seems paradoxical; one would expect that as time 
increases the distribution would get narrower, reflecting 
the fact that the observer is more confident about which 
eigenstate they actually have. This paradox is resolved 
when one realizes that the distance between the peaks 
increases at a rate greater than the increase in the width 
of the peaks. This time-dependent scaling of the distance 
between the peaks can be removed by changing variables 
to V = R/ {2y/2rjt). Under this transformation the prob- 
ability distribution is 



dVV[V) 




-(y-s)V(i/47t) 



(11) 



s=-J 



Now the FWHM of each peak scales as yJ\\Y2/{'^t) w 
0.83/V7t, which clearly illustrates that the probability 
distribution becomes sharply peaked about the D values 
for long-times. In Fig. [2] the dashed curve is a plot of 
Eq. ([11]) for £> = 5 and t = 47~^. All of the peaks are 
clearly distinguishable with a measured FWHM 0.418. 
This agrees with the prediction given by 0.83/y^. 

Returning to the non-symmetric sums, i.e. Eq. ([7l), for 
the remainder of this calculation, the purity is 



Tr 



The impurity, for a particular value oiV, is given by 

Mv,t) = (1- ^-°' ) (13) 



(14) 



The sohd curve in Fig. [5] is a plot of Eq. dH]) for D = 5 
and t = 4'y-^. To find the average value of impurity for 
a continuous measurement of Jz one must integrate this 
function over the scaled measurement record V, weighted 
by P{V) i: 



K{t,V)U{V)P{V)dV. (15) 



When using this result to calculate the measurement- 
induced evolution of the mean linear entropy one should 
use the kernel in Eq. (|14l) rather than the kernel in 
Eq. This is because it will reduce the numerical 

error that occurs when subtracting from one a number 
very close to one. 

To calculate this integral one must perform numerical 
integration. It would be useful to have a fully analytic 
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FIG. 2: The dashed curve is the probability distribution 
V{V), Eq. ([TTJ; the solid line is A(l/,f), Eq. shown for 

D = 5 at t = A-y'^. The dotted curve is K{V,t)M{V)P{V) 
multiplied by 20 (so that it can be viewed on the same scale). 
Also plotted are regions for the two eigenvalue expression for 
the impurity integral. 



expression for {L{t)) in the I?-dimensional case. Now 
we will make some approximations to obtain a simple 
analytic expression for Eq. ([T5|) . We will show that the 
approximations allow one to place a bound on the full 
expression. Then we will confirm our approximations 
with numerical simulations in Sec. IIII Bl and Sec. IIII Cl 

The first approximation involves truncating the state 
matrix to two eigenvalues and then renormalizing. We 
call this the two eigenvalue approximation. The moti- 
vation for the truncation stems from the following ob- 
servations. In the long-time limit {t ^ 7~^) the true 
probability distribution, V{V) — J^{V)P{V), is sharply 
peaked in D places; see Fig. [2] Also in this limit one 
finds that A(V", t), Eq. ([U]), is sharply peaked in {D - 1) 
places; the peaks are between peaks of V{V). Wiseman 
and Ralph have pointed out, for a qbit, that this arrange- 
ment leads to poorly purifying trajectories dominating 
{L{t)) 0. In a qbit the poorly purifying trajectories are 
those trajectories with eigenvalues of the same magni- 
tude {V{t) « 0) in the long-time limit. Physically this 
means that the measurement has not projected, or the 
filter has not decided, if the state is in 2 = ±i. In qdits 
there are {D — 1) regions like this, where V{t) « r for 
r £ [J — i, J — |, . . . , J -I- |, — J +5]. In any of these re- 
gions only two eigenvalues contribute significantly to the 
integrand. Consequently it is reasonable to truncate the 
state to the two eigenvalues in a particular region and 
renormalize. 



The effect of the two eigenvalue approximation is two- 
fold: it increases the purity of the state, and makes it 
possible to derive a fully analytic expression for the im- 
purity. The eigenvalues of truncated state matrix are 
(Ao, Ai)/(Ao + Ai)!!!], where Aq is the largest eigenvalue 
of p and Aa > Af, when a < b. It is easy to show that p 
is majorized [3^ by p2 (that is, p -< P2), meaning that 
the original state is more mixed than our two eigenvalue 
approximation to it. From the fact that the purity is 
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Schur-convex ^ it follows that Tr [p^] < Tr [pj] . This 
means the impurity of the two eigenvalue approximation 
is a lower bound on the true impurity: L{p) > L{p2). 
In the long-time limit it is reasonable to expect that 
L{p) « L[p2) from the arguments above. 

Now we calculate the impurity of the truncated state 
matrix {L2{t)). We split the integral in Eq. (|15p into 
regions such that the same two eigenvalues are the largest 
two eigenvalues for all values of V in each region. The 
integration is then performed in each of these regions and 
then the regions are summed. Figure ([2]) illustrates how 
the regions should be split for D — 5. After a coordinate 
transform (for example in region I it would be y = y — 
1/2) the integral in region one becomes 



9^1 



D V TT y_i cosh {i-ftV) 



-dV. 



(16) 



After a similar transformation the integral in region II 
becomes 



(17) 




^11 = 

The total integral is thus 

(L2(i)) = 25Hii + (i? - 3)5Hi. 



(18) 



For t ^ the integrands become sharply peaked about 
1/' = with negligible contributions to the integrals from 
the region outside the interval V d [— i], so that 5Hn — 
In this limit we have 



{L2{t))c^{D~l)--J — 




4'ytV'- 



The final approximation is arrived at by noting that 
in the long-time limit the distribution in the nu- 
merator is broad compared to the distribution in 
the numerator for t 3> 7^^- Thus we can make 
the approximation: ^J^dVe~'^'^^^ ' / cos\\{^f2rfV) ~ 

dV^/ cosh(-y27V^) — Tr/y^. The final expression 
for the impurity is thus 



(i2(0)LT = 



2{D-l) 



-It 



D 



(20) 



This is the analytic expression for {L{t)) that we set out 
to find, and is the key result of this section. 

Because previous work has used the exact qbit results 
we write it out for future reference: 



(i(i))qbit = ^J — 




-A-ytV'^ 



2 V TT 7„oo cosh(47iy) 
The long time limit of this expression is 



qbit- 



7re 



-7* 



LT 



dV. (21) 



(22) 



As could reasonably be expected for a qdit under the two- 
eigenvalue approximation (effectively a two- level system) , 
the resulting average impurity is proportional to the qbit 
impurity with the proportionality depending on the qdit 
dimensionality. In Sec. IIIIBI we examine the accuracy 
and validity of this approximation. 



B. Comparison with numerics 

In Fig. [3] (a) one can see that for t > 7"^ the impu- 
rity, found by numerically evaluating Eq. (jl5l) . is better 
approximated by Eq. (|T8)) than Eq. (|2ip . Moreover, for 
all times the qbit impurity Eq. ()21|) is a lower bound on 
the true impurity. However, when calculating the speed- 
up one is interested in times where ^ I, which are 
plotted in Fig.[3](b). 
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FIG. 3: A comparison of approximations to the evolution of 
{!/), for D = 3. (a) (starting from the upper- most curve): 
(L) i.e. Eq. H15|) shown by green circles; (L2) i.e. Eq. p8|) 
with the integrals evaluated numerically, shown by magenta 
upward pointing triangles; {I/)qbit i.e. Eq. (|2ip shown by 
solid black line, (b) (starting from the upper-most curve): 
(i2)LT i.e. Eq. (|20p shown by magenta downward pointing 
triangles; (L) i.e. Eq. (|15|l shown by green circles; {Z/)qbit-LT 
i.e. Eq. (f22ll shown by dashed black line; (I/)qbit i.e. Eq. ([2T|l 
shown by solid black line. Eq. (|18p is not plotted in (b) be- 
cause on this scale it is indistinguishable from (L). 



Eq. is not plotted in Fig. |3] (b) because on 

this scale it is practically indistinguishable from (L) 
at long-times (it approaches (L) from below). But 
{L2{t))LT, that is Eq. (pOj) . approaches (L) from above. 
Thus the approximations made to obtain Eq. (j20p from 
Eq. (|18p increase the impurity. This effect is eas- 
ily understood once one realizes that the approxima- 
tions amount to: throwing away parts of the integral 
Eq. (|19p : and over-estimating the integral by approximat- 
ing exp {—A'ytV''^) = 1. The next closest curve is that of 
the qbit long-time limit expression Eq. (HH), and below 
that is the qbit exact result Eq. (|21l) . 

From these curves we infer the following: the long- 
time expression for the two eigenvalue approximation, 
{L2{t))ijT, is closer to {L) than the qbit long-time 
limit; (I/2(t))LT is approaching (L) in the same way 
(L(t))qbit-LT is approaching (L(t))qbit- This is impor- 
tant, as previously in the literature it has been com- 



mon to use the qbit long-time limit expression as a lower 
bound for (L) any dimension [E, [l] . It is clear that for 
all times t ^ 7^^, {L2(t))irr is a better approximation 
to (L) than the qbit expression. Of course asymptoti- 
cally the expressions decay at the same rate. Therefore 
there is nothing in this work to suggest that the speed-up 
calculated in previous works are incorrect. 



C. Stochastic simulations 

We now compare the analytic solutions from Sec- 
tion |IlT3 and the exact numerics from Section fill Bl to 
stochastic non-linear trajectory simulations of Eq. ([l}. 
We used an Euler integration method with the following 
parameters: 5t = 1 x 10~^7~^; D = 5; ensemble size 
= 20. The ensemble size is small because we will plot all 
of the trajectories so that we may gain some qualitative 
understanding from them. 

Figure m (a) shows that the L of all of the trajectories 
is bounded from above by i for t 3> 7^^- Intuitively this 
effect can be understood from the two eigenvalue approx- 
imation; in the long-time limit a poorly purifying trajec- 
tory has two eigenvalues and thus these trajectories are 
upper bounded by L = ^. A more rigorous explanation 
is given in Sec. IIIIDI On this Figure we also plot the nu- 
merically calculated ensemble average impurity and the 
linear trajectories solution, Eq. (IT5|) . 

There is some difficulty in obtaining convergence for 
(L) in the stochastic simulations, even with very large 
ensemble sizes. The convergence problem is due to a 
small number of poorly purifying trajectories first noted 
for a qbit in Ref. 6]. These poorly purifying trajectories 
are evident in Fig. |4] (b) in that many trajectories touch 
the dot-dashed line aX L = i. For moderate ensemble 
sizes the whole region of possible impurities is filled by 
the ensemble even at i = 27^^. As noted above the 
major contributions to (L) are from the regions where 
V = r (see Fig. [2]); that is, for results V that are as far 
as possible from the most likely values of V [the peaks of 
'P{V)]. Thankfully one can use Eq. ([T5|) to calculate the 
evolution of the mean impurity. The advantages of this 
method are two fold: the computation time is greatly 
reduced, and it is exact, to numerical precision. 

Figure S] (b) shows that many trajectories also 
touch the dashed line corresponding to a scaling L ^ 
exp(— 47i). The reason for this clustering, and that at 
L = i, will be explained in Sec. IIIIDI 

D. The distribution of impurities 

In Ref. Jj6j Stockton, van Handel, and Mabuchi gave 
reasonable but non-rigorous arguments for a bound [their 
Eq. (41)] on the degree of spreading of trajectories (that 
is, on the width of the distribution of trajectories) for 
the model of Eq. ([TJ . We now use linear trajectories to 
rigorously derive the bounds that they found. In all cases 




Time (y-i) 



FIG. 4: Short time impurity dynamics for D = 5. We plot 
20 impurity trajectories (grey); the ensemble average of the 
trajectories (blue); the error bars for the ensemble average, 
±o"mcan, are smeared around the mean (light blue). Time is 
measured in units of 7~^. Fig. (a) Linear abscissa. The 
dashed and dot-dashed lines are approximate upper bounds 
on the L trajectories; see Sec. IIII Dl for further explanation. 
The solid black line is the linear trajectories solution for the 
qbit, that is Eq. (|21|) . The circles are (L) for D = 5, i.e. 
Eq. (|15|l . Fig. (b) the impurity is plotted on a logarithmic 
axis; this helps to see the approximately exponential decrease 
in impurity predicted by Eq. (|20|l . The approximate upper 
bound on the spread of L trajectories is the dot dash line. The 
dashed line corresponds to the most likely trajectories. Below 
the solid (magenta) line there is only l/D of the ensemble 
of trajectories at long-times. These lines are explained in 
Sec HTTD] 



it is easy to obtain an analytic expression, but the exact 
form of the expression depends on D. Therefore we quote 
only the long-time asymptotic scaling of the bounds. 

Central to our explanation of the bounds of the distri- 
bution of trajectories is the intuition gained by examining 
Fig. [51 where we have replotted portions of Fig. [5] on two 
separate axes. In Fig. [5] (a) 7^(V^, t) is plotted. This prob- 
ability density plot is significant in the current analysis 
because it shows which records, V , are likely (the peaks 
of V) and which are unlikely (the troughs). In Fig. [5] 
(b) the kernel A(y, t), that is Eq. ([14]), is plotted on a 
log scale. The key features of this plot are: 1) the peaks 
at V — — f,— ^j^jf which represent very impure tra- 
jectories and correspond to unlikely records, and 2) the 
troughs at F = — 1, 0, 1 coincide with some of the purest 
trajectories and correspond to minima of 7^(y, t). We will 
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elaborate on other relevant features of this graph below 
as each bound is explained. 
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FIG. 5: (a) probability distribution V(y,t), Eq. for 
D = 5 and t = 47~^. The shaded regions of the probabil- 
ity distribution contain 1/D of the total probability at long- 
times, (b) the solid blue hue is A{V,t), Eq. HH), plotted with 
logarithmic scale on the abscissa for D = 5 and t — 47^^. 
The red dashed line is a projection onto the abscissa of the 
minima of A(V^, t) a.t V — —1, 0, 1. The red solid line is a pro- 
jection onto the abscissa of the A{V, t) evaluated at V = ±2; 
notice that this projection is purer than the dashed line. The 
shaded regions in this figure depict the impurity of records 
that have |V| > J. As these regions concur with the regions 
in Fig. (a) they can be interposed as saying that only a pro- 
portion 1/D of the trajectories will have impurities smaller 
than those given by the projection A(±J, f). 



1. The upper bound: This is obtained by substituting 
V{t) = r (recall that r e [J-i, J-|,. . . ,J+^,-J+^]) 
into Eq. ([B]), normalizing, and then calculating the im- 
purity. This procedure gives the upper bound because 
these records correspond to the most impure trajecto- 
ries, that is peaks in Fig. [5](b). Because this procedure 
also corresponds to a a worst case scenario of the two 
eigenvalue approximation it can also be interpreted as 
the filter being unable to decided between the eigenval- 
ues A 1 and A 1 . For D = 5 the records of interest 

correspond to V = —-^i ^ (see Fig. [S] and Fig. [2]). 
In the long-time limit all solutions give 

Llt^I (23) 

which is also the bound predicted by Stockton and 
coworkers [l^. For small times there is some varia- 
tion between the bounds given hy V = r = ±^ and 
V — r = ±|. This is because V{V) is broad at short- 
times so that more eigenvalues contribute to the impurity 
when for example V = \. Conversely when V = ^ the 
is only one eigenvalue to the left of this record. Accord- 
ingly it is reasonable to expect that at short-times the 
bound given by ^ = ±i will produce a bound that is 
more impure than the bound given by V — ±|. This 
behaviour is confirmed in Fig. 2] (a), where the dashed 
line is I^ = ±1 and the dot dashed line is y = ±i. 



2. A pseudo lower hound: When V(t) = s for s G 
[J — 1, . . . , — J -I- 1], the record corresponds to the inner 
peaks of the probability distribution, as depicted in Fig. [5] 
(a) (the peaks at ±J are not included in this analysis). 
These peaks of 'P{V) are the most likely records and co- 
incide with some of the purest trajectories, the minima 
in Fig. [S] (b) . For D — 5 the peaks of interest are those 
corresponding to ^ = —1, 0, 1 in Fig. [SJ By substituting 
V{t) — s into Eq. ([6]), normalizing, and then calculat- 
ing the impurity, one obtains the lower bound without 
any approximations. These solutions can be obtained 
for any dimension D. By making a two-eigenvalue ap- 
proximation one can solve the evolution of the impurity 
analytically giving 

LLT--exp(-47t). (24) 

Equation (|24p corresponds to the lower bound quoted by 
Stockton and coworkers [16]. In their case, their bound 
appeared to be a true lower bound as their initial state 
was a coherent spin state (a collection of D/2 spin one 
half particles where D = 2 J + 1 with mean spin vector, 
J, aligned along the a;-axis in this case). This has a small 
population in |±J), and it is therefore unlikely in their 
case that \V\ > \J\. For our analysis we see clearly that 
the bound is not a true lower bound (see Fig. |4] (a)), 
because a maximally mixed state has equal populations 
in all eigenstates. However, because these records V{t) = 
j are the most likely records many of the trajectories 
touch or cluster around this bound. This will be further 
explored in point 4. below. 

3. Physically likely trajectories: In our case there is 
no lower bound, although a trajectory that purifies in- 
finitely fast is infinitely unlikely [59;]. It is possible to 
give a natural bound on the physically likely trajectories 
based on the probability distribution P{V). In the lower 
bound section above we did not include the peaks of the 
probability distribution found at = ±J. It is clear 
from Fig. [5] (b) that there is no corresponding minimum 
of the kernel A{V, t) in this region. When A(V, t) is eval- 
uated at V{t) — ±J one finds that it gives impurities 
smaller (see the solid red line in Fig. [S] (b)) than those 
which correspond to the inner peaks, [the dashed red line 
in Fig. [S] (b)]. From the shaded regions in Fig. [S] (a) it 
is apparent that only 1/D of the probability density is 
found for \V\ > J. This means that, in the long-time 
limit, only a proportion 1/D oi the trajectories will have 
purities greater than the bound we now present. Using 
the linear trajectory solution with V{t) = ± J and solv- 
ing for the impurity one obtains the bound. By making 
a two eigenvalue approximation one finds that the bound 
scales asymptotically as 

LLT~exp(-47i). (25) 

The scaling in Eq. (|25l) and Eq. (|24l) is precisely the same. 
The difference in the bounds is only apparent when the 
calculation is exact; see Fig. U] (b). The dashed line 
(black) is Eq. (1^^ . while the thin solid line (magenta) is 
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Eq. ([25]) (also see Fig. [6]). Obviously for a qbit (D = 2) 
the bound given by = ±^ corresponds to the median 
of the distribution of impurites, denoted by p{L, t)dL, at 
all times, and the mode at long-times. 

4. Distribution of Impurities: The distribution of im- 
purites can be calculated analytically for 13 = 2 6] . How- 
ever we will change variables to the log-impurity so that 
we may clearly see the featured discussed above. The 
required change of variables from Eq. (24) in Ref. [g'] is 
£ = logio L = logio [5(1 - ^^)] so that z = Vl - 2 x 10^. 
The distribution of log-impurites is thus 



t) = d£- 



e"^* In (10) cosh arctanh^(Vl - 2 x 10^ 



2y/Mjty/l - 2 X lO*^ 
X exp (-arctanh^(Vl - 2 x 10^)/47i).(26) 

In Fig. [S] this distribution is plotted at t = 27^^ (the 
red dashed curve). The peak of this distribution at 
£ — —0.301 becomes less prominent at long-times. This 
is when, as Wiseman and Ralph pointed out [6^], consid- 
ering the average log-impurity is a good way to find an 
approximation to the mean time to a fixed L. 

The poorly purifying trajectories affect (L) and a large 
number of these trajectories cluster around the upper 
bound (as seen in Fig. U]). This motivates the con- 
sideration of another measure of mixedness which de- 
emphasises the peaks of the distribution fsl, f6|. However, 
in this article, we will remain focused on (L) . 




-2 -3 -4 -5 

log-impurity 

FIG. 6: The distribution of log-impurities at t — 27"^. The 
blue solid curve is the distribution for D = 5. The red dashed 
curve is for D = 2, that is Eq. (|26p . The line with the circles 
is logiQ((Z/)) — —1.46 at t = 27"^, according to Eq. (fTS)) : the 
line with the squares is {£) — —2.41 at t — 27"^. Only a 
proportion 1/D of the trajectories have lower purity than (to 
the right of) the line with the crosses (plotted at £ = —3.1736) 
as explained in the main text. The two other features on the 
graph at £ — —0.3010 and £ = —2.87 correspond to the upper 
and the pseudo-lower bounds on the distribution of impurities 
described in the main text. 



For D > 2 the distribution of £ can be numerically 
calculated from the probability distribution ViV) and 
A(y, t) as follows. First one splits the distribution A{V, t) 



into regions {(—00, —j) 



-j, -j + 



. , [j, 00)}. These 



regions are chosen because £ is always monotonically in- 
creasing or decreasing in these regions; see Fig. [5] (b). 
This is required to perform the numerical interpolation 
required by the following procedure. In each region one 
discretizes the range of £ and then finds the values of V 
that correspond to this £. The probability distribution 
'P{V) is then integrated over this region to give d£p{£, t). 
This is then summed over all regions. 

In Fig. inithe solid curve is the probability distribution 
for £) = 5 at t = 27"^; this distribution was obtained 
using the above method. The upper bound and pseudo- 
lower bound described in the text above are quite ap- 
parent in distribution; they correspond to the two sharp 
features at ^ = —0.301 and £ — —2.87. It is also apparent 
that {£) is a more faithful central tendency measure than 
logj^Q(L) as the line corresponding to {£) is closer to the 
bulk of the distribution. 



E. Summary 

We conclude our study of information extraction 
from commuting measurements by summarising the key 
points. It is possible to characterise the average amount 
of information extracted using the average impurity. 
We found a simple expression for the average impurity, 
{L2{t))ijT ~ exp(— 7i), which is valid when t ^ 7"^. 
Additionally we found that in this asymptotic limit a 
portion {1 — 1/D) of the trajectories would be bounded 
by cexp(— 47^) < L < ^ for a constant c that we can 
determine for a particular D. It is the existence of rare 
but poorly purifying trajectories which explains the up- 
per bound of |, and which explains the difference be- 
tween the scaling of the mean and that of the lower 
bound, cexp(— 47^). Indeed, other measures of central 
tendency, such as the median, or the exponential of the 
mean-log, exhibit the same scaling as the lower bound, 
^ exp(— 47^), which is_also reflected in the mean time to 
attain a given purity [2] . 

Despite its limitations, we have used (L), which scales 
as exp(— 7^), as our measure of average information in 
the main text of the paper. Here we briefly comment 
on what would change if we were to have used a more 
faithful measure of central tendency such as the median. 
In Sec. IIVBI we found that the spread of trajectories 
for the QFT protocol was small for sufficiently frequent 
feedback. This feature is true in general for measure- 
ment in an unbiased basis The main consequence 
of this, for complementary continuous measurements, is 
the mean impurity is approximately equal to the median 
impurity. Thus we can directly compare the results ob- 
tained for (L) in the main part of the paper to the scaling 
for the median impurity, ^ exp(— 47^), for commuting 
measurement. The result, obviously, is a diminution in 
the speed-up offered by the latter by a factor of four: 
\{D + l)< S'^^'ii^" <\D'^- For smafl systems (D = 2 or 
3) this implies a slow-down, or at best a modest speed- 
up. But for large D the scaling with D implies that using 
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feedback to construct suitable measurements in an unbi- 
assed basis will beat the no-feedback protocol even in 
terms of the median impurity or mean time. 

IV. INFORMATION ACQUISITION IN 
DIMENSION D AND ITS RELATIONSHIP TO 
COMPLEMENTARITY 

A. Previous results 

In this subsection we briefly summarise the results of 
previous studies on information acquisition by measure- 
ment in a complementary basis. The original articles on 
the subject are found in Refs. 0, SBQ- The basic idea 
is to use quantum feedback control to continuously keep 
the eigenbasis of the state complementary to the mea- 
surement basis. By analogy with the intuition presented 
in Fig. [1] it is hoped that this procedure will enhance 
information extraction. 

It is easy to show, from Eq. ((l}, using Ito calculus 
that the information acquired [4I , as characterised by the 
change in impurity, about the state due to measurement 
of infinitesimal duration is ^ 

dL = -87rft{Tr [pXpX] ~ 2Tr [Xp] Tr [Xp^] 
+ Tr[pX]'Tr[p2]} 
-4V27dw;{Tr [Xp^] - Tr [pX] Tr [p^] }. (27) 

We wish to study the effect of measuring in a complemen- 
tary basis. In order to measure in a complementary basis 
through out the measurement process one must use quan- 
tum feedback control 13] to continually adjust the basis. 
We will not labour on those details here as they have been 
adequately discussed before 0, d, S [l3|- The required 
transformation of the measurement basis is X{t) X{t) 
such that — 1/\/D, for all x and i where \i) is an 

eigenstate of p and \x) is an eigenstate of X. Replacing 
all the X's by X in Eq. (|27|) considerably simplifies the 
expression to [1, Q 

dL = -87Tr [XpXp] dt, (28) 

or equivalently 

dL = -S-idt'^\XrA'^\r\c. (29) 

Observe that the unbiased condition which simplified 
Eq. p7)) to Eq. (|28)) leaves a permutational degree of 
freedom in Eqs. ([28|) and ([29|) because a permutation of 
the eigenstates of p does not affect = 1/Vd. This 

means we should optimise over permutations in the ex- 
pression dL — — 87Tr \X P^pPmX P}j^pPm\ dt, where Pm 
is a permutation matrix, to maximize the decrease in im- 
purity. To obtain a lower bound \dL\ of the optimal per- 
mutation the we average over all possible permutations: 
dL = -87S„Tr [XPlpPmXPlpPra] dt. The average 



gives a lower bound because the sum of a sequence is al- 
ways less than or equal to the greatest term in the sum. 

It was shown in Refs. [1, Q that a tight lower bound 
on Eq. (|28p . for the optimal permutation, is 

dL<~\t{D + l)L{t). (30) 

While the upper bound was recently found to be 

dL > -]ptD^L(t), (31) 

for D > 1 Equations ^ and 1^ hold for measure- 
ment in any complementary basis. Using these bounds 
one may compare the time it takes for a commuting mea- 
surement (i.e. Eq. (120^ ') to extract a certain amount of 
information (i.e. to attain (L) = e) to how long a com- 
plementary measurement may take to extract that same 
amount of information. That is we equate Eq. (j20p and 
the solutions of Eqs. (|30l) and (|3T|) . and solve for the ra- 
tio of icompiomontary and icommuto- Wc caU this ratio the 
speed-up in information acquisition. For tcommuto S> 7~^ 
the speed-up is bound by 

l{D + l)<S<^, (32) 

for L» > 1 0. 

While these results show the usefulness of feedback in 
rapid purification, there are some important issues that 
have not yet been addressed. First, does an explicit con- 
struction of such a complementary measurement strategy 
exist? Second, is actual performance of such a strategy 
0{D) or 0{D^)7 Third, what is the mechanism of the 
speed-up? The first and second questions are answered 
in Sec. IIVBI and Sec. IIVCI respectively. We defer the 
third question until Sec. |Vl 



B. Complementarity via quantum Fourier 
transform for D = 2, 3, 4 

In this subsection we explore implementing a contin- 
uous complementary measurement in systems of dimen- 
sion three and four. There are many unbiased bases, so 
to make our analysis concrete we must choose one. It 
has been noted elsewhere [l^ that for D > 3 not all un- 
biased bases are equally good at reducing the impurity. 
Nevertheless we choose a particular complementary, the 
one generated by the quantum Fourier transform (QFT) 
of the logical basis. The construction of the continu- 
ous complementary measurement then consists of: cal- 
culationally diagonalizing the state matrix; calculation- 
ally ordering (permuting) those eigenvalues to maximize 
|(c?L)| after a measurement; and applying the appropri- 
ate permutation and the QFT unitary to the state. This 
whole procedure is an example of quantum feedback con- 
trol ^] . The D-dimensional QFT can be represented by 
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the matrix 
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(33) 



where q — exp (2TTi/D). The matrix elements of T are 

T,, = -^exp(^^(rxc)^, (34) 

where r, c S [0, 1]. 

For D = 2 Jacobs' protocol [1] turns out to be equiva- 
lent to the QFT protocol. The transformation unitaries 
in both cases are 



exp(if Jj,) = -1= ( J J 



(35) 
(36) 



The effect of these transformations on dL is identical: 

^ = ^Jac-^^^lac = ^QFT-^-^QFT = BecauSe this 

scheme is equivalent to Jacobs we will not analyse this 
D = 2 case any further. 

For a system of dimension D = 3, the transformed 
measurement is 




(37) 



where q = exp(27ri/3). Taking the state to be p = 
diag(Ao, Ai, A2), where A^ > Ac when r < c and solv- 
ing the equation for dL, using the fact that L = 

Er#c,c#r ^rK, WC find 

dL = -S-idt^\XrA'^\r\a 

r.c 



-87dt ( -(A0A1+A0A2 + A1A2) ) (38) 



-8-fdt-Lit), 



(39) 



which coincides with the bound of Eq. (pOj) . The above 
expression can be easily integrated to give 



Lit) 



e 3 



^*L(0), 



(40) 



In Fig[7]the weights |Xr,cP between the ArAc's are plot- 
ted. The weights in this case are all equal, so it is obvi- 
ous that permuting the eigenvalues will not change the 
decrease of impurity. This explains why we did not have 
to find an optimal permutation. 

Now we compare the above calculation to numerics. 
The average impurity for the continuous complementary 




FIG. 7: The only possible permutation for D = 3. The solid 
lines between the eigenvalues represent the weighting of the 
multiplying factors in Eq. [38l That is, | = |Xr,c|^, Vr, c. 
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FIG. 8: Trajectory simulations of the impurity for a D = 3 
system subject to feedback every St = 10~^ 7"^ to keep the 
state and measurement eigenbases complementary. The en- 
semble size is 100. The dashed line is (L) for a continuous 
commuting measurement (D = 3) i.e. Eq. (|15|) . The solid 
line is the numerically calculated ensemble average of (L) for 
the complementary measurement; it is plotted in green; error 
bars are plotted at times (0.5, 1, 1.5, 2) x 7"^^. The analytic 
expression for the average impurity, Eq. (|40[) . is indistinguish- 
able from the results of the simulation. 



measurement protocol and continuous commuting mea- 
surement are calculated by stochastic simulations and 
Eq. respectively. Figure [5] shows that the stochastic 
simulations of the complementary measurement protocol 
give a deterministically decreasing (L) as predicted by 
Eq. (gni). 

We numerically investigated the effect of varying the 
frequency, l/6t, of the operation applied (feedback) to 
keep the state and measurement eigenbases complemen- 
tary in Fig. [3] The stochastic fluctuations in the simu- 
lations arise from the finite size of St in the simulations. 
This effect can also be found in simulations of the Jacobs' 
qbit feedback protocol '6^. As (5t — > the numerically 
calculated average impurity for complementary measure- 
ment approaches Eq. (HO]) . For 5t ^ 10~'^7~^ numer- 
ical simulations are indistinguishable from the analyti- 
cally calculated (then interpolated) speed-up, the solid 
line in Fig. [3] 

We now consider complementary measurement using 
the QFT for the case when D = A. The increment for 
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FIG. 9: Short time numerically calculated speed-up for D — 
3. The dashed line is the analytically calculated asymptotic 
speed-up S = 8/3 — 2.6. The solid (magenta) line is the ideal 
speed-up; calculated numerically from equations Eq. H15I) and 
Eq. (|40p . The circles and triangles are numerical calculations 
of the speed-up with finite St — 10~"^7~^ and 10~^7~^ respec- 
tively; the ensemble size is 100. As iSt — > the numerically 
calculated speed-up approaches the solid line and the error- 
bars vanish. 



the impurity is 

r.c 

= -8jdt (AqAi + iAoAa + iAiAa 

-fAoAs-f-AsAg + AiAs). (41) 

Unfortunately, it is not possible to factor out an ex- 
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(a) (b) (c) 

FIG. 10: (a) The weights \Xr,c\'^ in Eq. (jil]) after diagonaliza- 
tion and descending ordering, (b) The optimal permutation 
for D = 4. (c) The worst permutation for D = 4. The 
solid line represents maximum weighting of multiplication in 
Eq. (|41l) . the dashed lines represent minimum weights corre- 
sponding to i. 



pression for the impurity on the RHS of the above equa- 
tion. It is obvious from Eq. (|4ip that permuting the 
eigenvalues will affect \dL\. Our task is to maximize 
the decrease in dL. From Eq. (|^1]) it is possible to in- 
tuit the form of the best and worst permutations, which 
are depicted in Fig. [10] An optimal permutation is 
Popt = diag(Ao, Ai, A3, A2). This permutation is not 
uniquely optimal; the permutation diag(Ai, Aq, A2, A3) is 
also optimal. One of the worst permutations is Pworst — 



diag(Ao, A3, Ai, A2). Unfortunately knowing the optimal 
permutation does not help us simplify the expression. 
Because we cannot solve this case analytically, we invoke 
the procedure developed in Refs. to find the bounds 
on \ {dL)\ (and hence S) for the QFT protocol for all D. 
These bounds are calculated in Section HV CI 

Figure [TT] depicts the numerically calculated speed-up 
for D = 4. The speed-up is larger than that of the D = 3 
case [6l| . Furthermore, for reasonably frequent feedback 
the numerically calculated speed-up lies between the an- 
alytical bounds predicted by Eq. ([50]) and Eq. ([5T|) . 
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FIG. 11: Short time numerically calculated speed-up for D = 
4. The dashed lines are the bounds on the speed-up given 
by Eq. (|44p and Eq. (|48p . The circles, squares and triangles 
are numerical calculations of the speed-up with finite St — 
10~^7~^, 10~^7~^, 10~^7~^ respectively. To keep the error 
bars roughly the same size the ensemble size was varied. The 
ensemble sizes for circles, squares and triangles are 100; 400; 
and 6, 400 respectively. 



We also briefly compare the complementary measure- 
ment protocol to the rapid measurement (RM) algorithm 
of Ref. in Fig. [T^l The decrease in the ensemble av- 
erage of the impurity seems to be of the same order. 
However, the trajectories in the RM case have a large 
variance. This is because the trajectories are not difFer- 
entiable. By comparison, even for a finite dt the L in 
complementary measurement is very close to determinis- 
tic. Although L is still stochastic, it is also diflterentiable 
for St — >■ dt so the noise in the simulations is reduced. 



C. Bounds on the QFT complementary 
measurement protocol for all D 

1. Lower hound 

From Eq. ([5^ it is clear that the lower bound on the 
optimal speed-up for measurement any complementary 
basis is 5 = |(£) -I- 1). Nevertheless it is interesting to 
work this out explicitly for a particular complementary 
basis. The method we use to find the lower bound was 
first presented in Ref. This method uses a fictitious 
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FIG. 12: A comparison between the complementary mea- 
surement protocol and rapid measurement (RM) feedback (of 
Ref. [3) for L» = 4; = 10"*7"^ and an ensemble of 20 
trajectories in each case. The dashed lines are the analyt- 
ically calculated bounds on dL for complementary feedback 
[Eq. (|44)) and Eq. (|47)) ]. The dot-dashed line is the ensemble 
average for L with QFT feedback. The solid (dark green) line 
is the ensemble average for the RM feedback. The light brown 
lines are the trajectories for the RM feedback. The analytical 
bounds suggest that both feedback strategies have the same 
order of speed-up. Inset figure: all of the complementary 
measurement trajectories (grey lines) lie within the predicted 
bounds. The trajectories are very nearly deterministic. 



state pf with the same purity as the state under con- 
sideration, i.e. L[p] — L[pf]. The fictitious state pp is 
defined as pp = diag(l — A,S, ..,S) where 1 — A is the 
largest eigenvalue of pp and S = A/{D—1). The equation 
for dL is invariant under permutations of the eigenvalues 
Pf- Substituting pp into the RHS of Eq. (P^t gives 



2(1 -A) 



A 



(D-l) 



E 



This expression is further simplified by explicitly calcu- 
lating the matrix elements. For the QFT we have 



D-l 



1 



4^ ^ 4sin2 (Trm/D) 



(i^2_i)/i2, (42) 



using the identity YJLA csc^i-Km/D) = (L>2 - l)/3 [3£ 
Also, 



E 



D + l 
12 



{D-l){D-2) 



(43) 



using the above identity and the fact that csc^(7r(_D — 
1)/!?) = csc^(7r/D). After some simplification the final 
expression is an upper bound on the decrease in impurity: 



dLp < --{D + l)jdtL{t), 



(44) 



which implies that the lower bound on the asymptotic 
speed-up is 



SF>-iD + l), 



(45) 



as found in Refs. [1, Q for a general complementary mea- 
surement. 



2. Upper bound 

To find the upper bound on the speed-up for the QFT 
protocol we again uses the method of Ref. [9]. This 
instance requires a different fictitious state p2, again 
with L[p] = i[p2], which was termed the binary dis- 
tribution. The binary distribution is defined as p2 — 
diag(l — A', A', 0, 0) where 1 — A' is the largest eigen- 
value. It is known that dL is most sensitive to permuta- 
tions of the eigenvalues of y02 @- 

Substituting p2 in to Eq. (|29l) gives only two terms. 
Now there is a choice about where to situate the two 
eigenvalues so that dL is maximized. By calculating the 
element |A"r.cP, using the QFT, we find that 



2 1 



cos 



D ' 



(46) 



This allows us to optimisze over the permutations by find- 
ing max \P^Xr,cPm\'^ ■ The largest element of the matrix 



IX, 



on 



1/2 (1 



D , 



From this it is clear that 
dL will be maximized provided the two eigenvalues from 
P2 are in succession. The change in the impurity for the 
binary distribution is 



dLo = 



-8-/dt{2XaXi) 
271^ 



-A-fdt 



1 



Thus the asymptotic speed-up upper bound is 

4 

5*2 < — 



L2(t). (47) 



For I? ^ 1 one finds 



^2 = 



2D' 



2tt 



Q.2D' 



(48) 



(49) 



In Fig. [13] we depict the conjectured optimal and 
worst permutations of eigenvalues of a D dimen- 
sional state p. The optimal permutation is p = 
diag(Ao, Ai, A3, As, • • • ,A6,A4,A2). 

In Fig.[T4]the numerically calculated asymptotic speed- 
up as a function of D is plotted. As expected, it is within 
the lower bound of Eq. (I45|) and the upper bound of 
Eq. (|48p . Further it confirms a very nearly quadratic 
speed-up: the fit shown is S" = 0.189152-]- 0.109i:' -I- 0.248. 
If, however, we choose to only fit the quadratic term, the 
fit is 



S = 0.191)2 



(50) 



For D > \Q this is a good approximation and is very close 
to Eq. (|49p . which was derived for long-times and large D. 
The numerically calculated asymptotic speed-up is only 
a multiplicative constant away from the ultimate upper 
bound on all complementary basis purification protocols, 
the dashed line {S = 0.5D'^), i.e. Eq. 
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FIG. 13: a) Weighting factors in Eq. after the feedback 
has diagonalized and ordered the eigenvalues in p. b) The 
conjectured general optimal permutation, c) The conjectured 
worst permutation. The weights are only shown from the Ao 
perspective. 
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FIG. 14: The asymptotic speed-up for reaching a given level 
of purity as a function of the system dimension. The solid line 
the fit in Eq. ([50|) . The dashed lines are the upper, Eq. (|48|l . 
and lower, Eq. (|45|) . bounds derived in the text for the QFT 
protocol. The dot-dashed line is the upper bound on all unbi- 
ased basis purification strategies, that is the RHS of Eq. p2)) . 



D. Discussion 

The results presented in sections Sec. IIV Al -Sec. IIV CI 
raise a number of questions which we now elucidate fur- 
ther. 

Because of symmetries, for dimensions two and three 
permuting the basis of p was unnecessary to maximize 
dL. The intuition is, for maximizing information gain 
the complementary measurement marginalises the signif- 
icance of which eigenvalue is the largest. For Z3 > 3 we 
must permute the basis of p, after diagonalization and 
before applying Tqft, to maximize dL. The permuta- 
tions applied do not increase the signal to noise ratio of 
the measurement unlike the rapid measurement proto- 
col of Ref. [5]. In Sec. |IVB] and Sec. IIV C 21 we argued 
that the permutations should be chosen to maximize the 
product of the two largest eigenvalues in Eq. (1^^ . While 
this recipe is operationally sound it lacks physical in- 
sight. For example, this recipe does not explain what 
the mechanism of the purification is and why the permu- 
tations can dramatically affect the rate of purification. 
These questions are addressed in Sec. |Vl 



The final open question concerns the efficacy of differ- 
ent complementary bases to change the rate of purifica- 
tion. In previous work Lund and Jacobs suggested that 
not all complementary bases were equally good at en- 
tropy reduction ^9|]. That is, they claimed that one par- 
ticular unbiased basis can reduce entropy more quickly 
than another. If this were true then the S = 0.2D^ upper 
bound we derived from the QFT basis would not neces- 
sarily be the true upper bound, which could be closer to 
S — Q.5D^ We numerically investigate this question in 
Sec. Ell 

V. WHY COMPLEMENTARITY IS NOT 
ENOUGH 

How is it that rotating to an unbiased basis provides 
any speed-up? Naively following the reasoning applied to 
qbits in Refs. [2-4] suggests that making each outcome 
equally likely {i.e. maximizing the surprisal [l|) one max- 
imizes the average amount of information a measurement 
extracts. From this one might argue that the speed-up 
observed can be explained by the Z3-dimensional version 
of this argument. However, it is not clear how this ar- 
gument would explain why the arrangement of eigenval- 
ues in p (permutations) are important to attain the best 
speed-up. To address this, we examine a phase space rep- 
resentation of the optimal and worst permutations that 
are schematically depicted in Fig. [T3] 

Defining a phase space picture for a discrete variable is 
not trivial. In Ref. [36| it was shown that a spin Wigner 
function W{9,(j)) can be defined in terms of Clebsch- 
Gordan coefficients and spherical harmonic functions. 
This spin Wigner function is a pseudo-probability dis- 
tribution on the Bloch sphere, with 9 and (f> the usual 
Euler angles. The spin Wigner function is a little counter- 
intuitive; for example, unlike the original Wigner func- 
tion, W{x,p), for position and momentum the marginal 
distribution for (f> is not the true phase distribution P{4>). 
However, for large D the marginals are a good approx- 
imation to the true phase distribution [37i] ■ We plot 
the Wigner function using the equal-area projection (de- 
scribed by co-ordinates <j) and J cos 9). 

The conjectured best and worst permutations were de- 
picted on a ring in Fig. 1131 explained in terms of the 
weighting factors jX^^cP- We now have an intuitive un- 
derstanding of the angle around this ring as being the 
phase (f) in the angular representation of spin states. The 
eigenstates of p in the measurement (Jz) basis are \r) 
(Dicke states [3^|). These are transformed by the QFT 
to the states \4>r), where (pr :— ^(J — r). The states 

\(j)r) = '^^m=-j^^Pi~^''^'f'r)\J,'m) are equivalent to 
the Pegg-Barnett phase states [H, . If the QFT were 
an easy operation in some physical system then our pro- 
tocol would be a procedure for rapidly preparing a Pegg- 
Barnett phase state. We note that in some spin systems 
it should be possible to construct the desired unitary ^ij . 
Consider the long-time-limit state. The worst case, for 




FIG. 15: An equal area plot of the Wigner function of an 
equal mixture of two eigenstates after the QFT, D = 10. a) 
the worst permutation cf. Fig. [13] (c). The dashed vertical 
lines correspond to the locations of the peaks of the D trans- 
formed eigenstates. The dotted box around the peak at <j!> = 
corresponds to the bulk of the Wigner function for this state, 
b) The conjectured general optimal permutation, cf. Fig. [13] 
(b). 



purification purposes, is when the eigenvalues are equal. 
Under this two eigenvalue approximation the worst per- 
mutation corresponds to putting the second largest eigen- 
value, Ai, the maximal distance away from Aq in phase 
space (0 — ±7r). The spin Wigner function for this con- 
figuration of the mixture is plotted in Fig. [15] (a) for 
D = 10. The optimal permutation is when the two 
largest eigenvalues are next to each other in phase space; 
see Fig. [13(b). It is now apparent that the schematic dia- 
grams in Fig. [12] represent a slice through the unwrapped 
Bloch sphere and the positions of the eigenvalues are their 
arrangements of the phase states corresponding to the 
original eigenvalues. This phase space picture suggests 
it might be possible to explain why these are the best 
and worst cases and how the speed-up is generated. To 
explain these two features we move away from exact cal- 
culation of Wigner functions and move to a schematic 
representation of the Wigner function. 

In Fig. [TH] (a) we represent schematically the bulk of 
the Wigner function by a rectangle of width 2tt/D (cor- 
responding to the (j) coordinate) and height D — 1 — 2J 
(which corresponds to the J cos 9 coordinate). In what 
follows one may loosely think of the rectangle as rep- 
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FIG. 16: Schematic Wigner functions for mixtures of two 
phase eigenstates for D = 20. (a) The worst permutation 
cf. Fig. [13] (c); we have rotated the entire sphere by n/2 
about z to simplify the explanation of the speed-up. b) The 
worst permutation after a positive measurement result (mod- 
ulo it/2). The total area is constant in figures (a) and (b), 
and is equal to Att. (c) The optimal permutation cf. Fig. [13] 
(b). (d) The optimal permutation after a positive measure- 
ment result. The total area is reduced in figure (d) to Stt, 
which gives rise to a purifying effect (the purity is inversely 
proportional the area underneath the Wigner function). The 
striped region (f) £ (— 7r/_D,0] denotes that the red and green 
rectangles are overlapping in this region. From these figures 
it is apparent that the circles in Fig. [13] can be thought of as 
showing (/!>. 



resenting the uncertainties in an observer's knowledge 
about two conjugate variables; according to the Heisen- 
berg uncertainty relation (HUR) the area of the rectangle 
must be constant for pure states. Now consider the effect 
of a weak measurement of J^. In our protocol this is of 
infinitesimal duration, but here we exaggerate the effect 
to illustrate our point. In Fig. [16] (b) we have taken the 
result of the measurement to be positive. Because the 
result was positive, the Wigner function does not have 
much support on the lower part of the plot. As the pos- 
itive result contained information about the distribution 
of j, the uncertainty in this variable is reduced. In keep- 
ing with the HUR, the conjugate variable {(/>) suffers an 
increase in variance. For Wigner functions the purity is 
proportional to the inverse of area under the function 
[42]. Here there has been no change in the total area 
and hence no change in the total purity. Because this is 
only a heuristic for understanding protocol the previous 
statement is not entirely true; the final paragraph of this 
section will explain the actual result. 

In Fig. [16] (c) the optimal permutation in the QFT ba- 
sis is depicted. The total area of the two rectangles before 
the measurement is ~ An. After a positive measurement 
result the phase distributions for the two eigenvalues sig- 
nificantly overlaps. The total area is now ~ Stt; this re- 
duction in area leads to an increase in purity; see Fig. 1161 
(d). Now the intuitive understanding for the speed-up 
and the permutation sensitivity is apparent. The per- 
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mutations are important so that the large eigenvalues in 
p may bleed into each other after a measurement. The 
bleeding is due to measurement backaction in the variable 
conjugate to J. It is the reduction in area this bleeding 
effects that causes the purification; but it only works if 
the largest two eigenvalues are adjacent in phase space. 
Although this picture is crude it captures the essence of 
the protocol. 

So far we have given an intuitive explanation of the 
mechanism of the purification. Now we provide an in- 
tuitive explanation for the speed-up. Consider the inte- 
grated measurement result 



/ dR{t) = / ^X{s)ds + [ 
Jo Jo Jo 



dW{s) 



(51) 



this expression contains a term corresponding to the sig- 
nal of interest (the first term) and a term representing 
the noise (the second term). We may define the signal to 
noise ratio as SNR = signal^ /noise^. In this case we have 
signal^ = AjAt^{X^) and noise^ = {{AW{t))^) = At. 
Thus the ratio becomes SNR = 47 At {X'^y Now we wish 
to estimate the time taken to evolve from the Wigner 
functions in the first row of Fig.[Tn]to those in the second 
row. Note that an observer gains one bit of information 
about Jz between the first and second rows. From the 
Shannon-Hartley theorem, bits = log2(l -|- S/N), we in- 
fer that this implies SNR = 1. Now we solve for time 
and take X = consequently (j^) (_D^ — 1)/12 (see 
appendix C in Ref. 0) to find 



At 



(52) 



Recall that the purity is approximately inversely propor- 
tional to the phase space area. The initial purity in 
Fig. [in] (c) is Pi oc l/47r and we know the state is an 
equal mixture of two eigenstates which means the purity 
equals one half. This allows us to determine the propor- 
tionality constant to be 27r, so = ^. The final purity 
is thus Pf = 2/3 which means the change in purity is 
AP = 1/6. Consequently 



lAF^i^i/e^^i 

P7At 2 3/(L>2_l) 9^ 
This can be compared to Eq. (H^ that is 



1 dL 
Ljdt 



2D^ 



5^^ 



(53) 



(54) 



Considering the crudeness of the arguments we have em- 
ployed, there is surprisingly good agreement between 
Eq. (|53p and Eq. ([54]) . This gives weight to the intuition 
that the effect of purification for the QFT protocol comes 
from measurement back-action and the role of the per- 
mutations is to maximize the back-action by placing the 
eigenvalues close in phase space. 

We now return to the worst permutations for purifi- 
cation. All one can conclude from the equi-likelhood 



of all outcomes is that there should be a speed-up of 
at least two (in a discrete outcome measurement model, 
which is equivalent to our current measurement model). 
Let us confirm this intuition now. Returning to the two 
eigenvalue approximation, consider the worst permuta- 
tion. It is possible to determine the upper bound on 
speed-up resulting from this permutation by substitut- 
ing r- c = D/2 into Eq. (gH]). We find that S = 2, which 
confirms the intuition above. This is why complemen- 
tarity is not enough to guarantee maximal information 
extraction. 



VI. ARE ALL COMPLEMENTARY BASES 
EQUALLY GOOD FOR ENTROPY REDUCTION? 

We now address the question raised by Jacobs and 
Lund |19i: are all unbiased bases equally good for rapid 
purification? Their answer was no. However at the time 
of their analysis the role of permutations was not clearly 
understood. Thus it is worthwhile to re-examine this 
question. Due to rotational symmetry of the unbiased 
bases and the permutational symmetry of density opera- 
tor the answer in the cases D — 2 and 3 is yes. However 
for _D = 4 it is easy to find a counter-example to this 
trend. For example measuring in any of the four mutu- 
ally unbiased bases (MUBS) ,43] gives 

dL = -8-fdt'^\Xr,c\'^XrXc 

r,c 

= ~8-/dt [2A0A1 + i(AiA2 + A0A3) + 2A2A3 

+ 0(AoA2 + AiA3)], (55) 

here X = AljjzMi where Mi is a transformation to one 
of the MUBS. (Explicit expressions of the five MUBS 
when D = A are given in Ref. for example. For 

convenience we have reproduced them in a footnote [6^.) 
As before, we may use p2 to obtain a lower bound on dL 
(and hence an upper bound on the speed-up). Doing so 
gives dL = —8^dtL{t) which implies a speed-up of S* = 8. 
This saturates the upper bound on the speed-up found in 
Eq. (|32p . The factor of two improvement over the QFT 
can be understood from the difference in weights of the 
largest terms AgAi between Eq. (|4T|) and Eq. (f55|) . 

To get more intuition about the purification process 
using the MUBs transform we look at the Wigner func- 
tion of the four states in one of the four MUBs, which for 
convenience we denote {|0), |1), |2), |3)}, in Fig. [HI The 
states plotted in Fig.[T7](a) and (d) look like the states in 
the QFT basis (Pegg-Barnett phase states). The Wigner 
functions plotted in Fig.[T7](b) and (c) are however quite 
different — they contain a hole which is more negative 
than the dips in Fig. [17] (a) and (d). The states plotted in 
Fig. [T7] (b) and (c) can be said to be highly non classical 
because of this. 

We now compare the MUB transform to the QFT us- 
ing the Wigner function, in an attempt to get an intu- 
itive explanation of the advantage of the MUB transform. 
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FIG. 17: An equal area plots of the Wigner function resulting 
from transforming four different states to the Mi unbiased 
basis {D = 4). (a) the transform of po = <liag(l, 0, 0, 0) (b) 
po = diag(0, 1,0, 0) (c) po = diag(0, 0, 1, 0) and (d) po = 
diag(0, 0, 0, 1). Unlike the unbiased basis generated by the 
QFT not all of the resulting states are like the Pegg-Barnett 
phase states; see for example figures (b) and (c). 



To make the analysis simple we restrict to states which 
have two large eigenvalues only (as we did above). From 
Eq. (|55l) we see that permuting these states in the logical 
basis so that they correspond to |0) and |1) (or |2) and 
1 3)) is the optimal thing to do. Recall that permuting the 
two largest eigenvalues to |0) and |1) was also optimal for 
the QFT protocol; see Fig. [18] (a). There we found that 
doing so created two peaks close in phase space, so the 
purification could be explained by measurement backac- 
tion. In Fig. [18] (c) the Wigner function of the optimally 
permuted state for the MUB transform is depicted [it is 
an equal mixture of Fig. [T7] (a)and (b)]. Here we can 
no longer attribute the purification mechanism to back- 
action. In fact it seems as though the purification effect 
comes from distinguishing the two peaks in the Jz distri- 
bution. The worst QFT permutation also corresponds to 
a poor MUB transform as evidenced in Fig. [18] (b) and 
(d). 

Interestingly there is a permutation that is worse for 
the MUB transform than the permutation depicted in 
Fig. [18] (d). It is p = diag(l,0, l,0)/2; the transformed 
state of this permutation is depicted in Fig. [TH](e). This 
permutation results in no purification at all! We can see 
this by substituting Aq = ^ and A2 = ^ into Eq. ([55]) to 
get dL = 0. This complete lack of purification may have 
application in state stabilisation of states with arbitrary 
purity [H]. 

We would like to know if the upperbound on the speed- 
up is saturated in all dimensions. To answer this question 
we resort to a numerical search for D G [2, 10]. We do 
not claim that our search is exhaustive. Recall that the 
stochastic simulations show that speed-up predicted by 
P2 is close to the achievable amount. It seems reasonable 
to assume this to be true in other complementary bases. 
This greatly simplifies the analysis by making the opti- 
misation of permutations superfluous. Because of this 
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FIG. 18: A Wigner function comparison between the QFT 
(top row) and MUB transform (bottom row) for D = 4. 
We plot an equal mixture of two eigenstates with the op- 
timal permutation (first column) and a poor permutation 
(second column) for each transform. In figures (a)-(d) the 
black dashed line is the zero contour, (a) the optimally 
permuted state for the QFT - p = diag(l, 1, 0, 0)/2; (b) 
the worst permutation for the QFT - p — diag(l, 0, 1, 0)/2; 
(c) the optimally permuted state for the MUB transform - 
p — diag(l, 1, 0, 0)/2; (d) a poor permutation for the MUB 
transform - p = diag(l, 0, 0, l)/2. Note that cases (b) and (d) 
are actually identical. Unlike the QFT case, the mixture of 
optimally permuted MUB states has two peaks which could 
be distinguished by a measurement of Jz. (e) The worst per- 
mutation i.e. p = diag(l, 0, 1, 0)/2 after the transform to the 
Ml unbiased basis, which results in no purification. 



we may numerically search for over all unbiased bases for 
the basis which has the largest element |Xr,cp. Convert- 
ing the element |^r,cp to a speed-up gives the following 
trend. The speed-up for even D attains the upperbound, 
i.e. -D^/2, and the speed-up for odd D equals (D — 1)^/2, 
as seen in Fig. [19] In either case this is much larger than 
the speed-up of the QFT protocol which was S = 0.2_D^. 
For -D ^ 1 it is reasonable to believe that the achievable 
speed-up does indeed scale like the predicted S — 0.5D^ 



VII. INFORMATION ACQUISITION IN A 
REGISTER OF QBITS AND ITS RELATIONSHIP 
TO COMPLEMENTARITY 

Finally, we consider complementary measurements of 
a register of n qbits, where each qbit is independently 
and weakly (or continuously) measured, as introduced in 
Ref. 0. Instead of one observable X, there are now n, 
given by X^'') = /(i) (g) ® . . . cr^'') . . . (g) /("), where r 
labels the rth qbit. The SME describing such a measure- 
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FIG. 19: The asymptotic speed-up for reaching a given level 
of purity as a function of the system dimension. The dashed 
lines are the upper bound and lower bounds on all unbiased 
basis purification strategies, i.e. Eq. ((32]). The dotted line is 
the maximum achievable speed-up for the QFT protocol i.e. 
Eq. H48|) . The solid circles are the results of a numerical search 
over unbiased bases to maximize the asympototic speed-up. 
See the text for more details. 



ment is 
dp = ^2ndtV 



X 



(r) 



p. (56) 



Here we define the complementary observable to be 
X^'''") = PmTX^'''^T'^Pl^. As before, the T's are condi- 
tional unitaries that introduce the unbiasedness (between 
p and X^''^), and the PmS are the permutation operators. 

It was found that the upper bound on the impurity 
was 



L^^l{t) = exp{-8Knt/[D-l])L{0), 



while the lower bound was 

-(n) 



L\^'^{t) = exp(-8Knt)i(0). 



(60) 



(61) 



From these equations we can infer the following bounds 
on the asymptotic speed-up factor 



2n 
D - 1 



< S <2n. 



(62) 



For all > 3 the lower bound on the speed-up becomes 
less than unity. For large n the slow-down implicit in the 
lower bound is ^ n2~"+-'^. (The interested reader should 
also see the related study of Hill and Ralph (ilj.) 



Complementarity via the quantum Fourier 
transform 



The combined state of the n qbits exists in a D = 2"- 
dimensional Hilbert space. 



A. Commuting measurements analysis 

Here we will not analyse the commuting measurements 
of a register of qbits with the same detail as we did for 
the qdit in Sec. IIIIl Instead we rely on a result from 
Ref. Q where an expression for the evolution of the av- 
erage impurity of a register by undergoing a continuous 
commuting measurement was found. In the long time 
(LT) limit the impurity is 



(57) 



We will use this expression in this section to calculate the 
speed-up, so only the asymptotic scaling is important: 



B. Complementary measurements 

It was shown in Refs. ^ that the change in impurity 
for a register of qbits monitored in an unbiased basis is 



n 

dL = -8Kdt^Tr [a(''V^^''^ 



(58) 



(D-l) 



-S^dtY^ ^ l^grA.A,. (59) 

r=l i,j=0 



From Eq. (I62p it is not clear if measurements in a com- 
plementary basis provide any benefit in information ex- 
traction. In this subsection we present some progress 
towards answering this question. 

We begin our analysis by considering the speed-up 
when n = 2. The bounds on the speed-up are 4/3 < 
S < A. In order to see if the upper bound is achievable 
we perform stochastic simulations. In Fig. [20] we nu- 
merically determine the advantage of the QFT feedback 
protocol in a register of two qbits. The permutation used 
for simulating the evolution was p = (Aq, A4, A3, Ai). 

In Fig. [20] (a) the QFT protocol clearly does not sat- 
urate the bound on (L) given in Eq. (|61l) . One may 
randomly permute the eigenvalue arrangement at times 
6t before applying the QFT feedback to determistically 
achieve the lower bound of Eq. (|62t . From Fig. [20] (b) it 
seems as though the speed-up is asymptoting towards 2 
rather than the upper bound specified by Eq. (|62]) . i.e. 4. 
For comparison, we note that the speed-up found for the 
locally optimal rapid measurement protocol in Ref. Q 
was S'rm ~ 1-4 in the long-time limit. 

The structure of unbiased bases for a register of qbits 
is quite complicated fiB - lirj . Even for the QFT the op- 
timal permutation is not obvious for n > 2. It is for 
this reason we plot the two eigenvalue approximation to 
the asymptotic speed-up for the QFT in Fig. [51] as a 
function of n. The values were obtained by finding the 
largest term in Eq. (|59p . This term will be denoted by 

ATmax ■— maxij- J2r "^^"^ feedback places the two 

eigenvalues at imax and jmax- Thus X^ax is proportional 
to the speed-up. Curiously the speed-up, 

5 « 2, (63) 
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FIG. 21: Bounds on the asymptotic speed-up for a register 
of qbits. The dashed lines are the upper and lower bounds 
derived in the main text Eq. (|62p . The solid line is the largest 
element of Xmax = maxi,j J2r ^he element Xmax is 

conjectured to be a good indicator for maximum achievable 
speed-up with the QFT in a register of qbits. 



FIG. 20: (a) Impurity for a two qbit register: dotted line 
- a continuous measurement; solid green line - average im- 
purity for feedback; solid grey lines - individual trajecto- 
ries grey lines; dashed magenta lines - upper [Eq. (|60|) ] and 
lower [Eq. (|6ip] bounds, (b) Speed-up solid line; lower bound 
dashed line; upper bound not shown as it well above the nu- 
merical result. The ensemble size is 20 and feedback was 
applied at times St = 1 x 10^'^7~^. 



appears to be independent of the size of the register. It 
is not clear if this is true for any unbiased basis in a 
register. Our arguments from Sec. IVII indicate it might 
be a multiplicative constant higher. 



VIII. DISCUSSION 

Prior to this work it has been shown that it is possible 
to speed-up the extraction of information from a quan- 
tum system using a continuous complementary measure- 
ment [1, H, 1^ . In this paper we have given an explicit 
method for constructing such a method, using the QFT, 
that achieves a speed-up of 5* = 0.2D^ for a qdit. This 
explicit construction allowed us to show that comple- 
mentary between the state and measurement observable 
is not enough to guarantee maximal information extrac- 
tion. Choosing the right permutation of eigenstates, be- 
fore Fourier transforming the state - see Fig.[T31 is needed 
to guarantee maximal information extraction, which is 
contrary to what one might expect from classical infor- 
mation theory [l|. The effect of the permutation can 
change the speed-up in information extraction from its 
minimal value 2 to 0.2D^. We also argued that in the 



case of the qdit, the maximum possible speed-up, pre- 
dicted in Ref. Q, 5 ~ should be achievable. One 
interesting upshot of our investigation (see Sec. IVI|) is 
an explicit example of using an unbiased basis for state 
stabilisation jl9i |. 

In this paper we demonstrated that measurement back- 
action is fundamental to the purification process for the 
QFT protocol. We did not, unfortunately, find an infor- 
mation theoretic explanation of the measurement process 
as was found in Ref. [i^. It is not clear, at present, that 
the bounds presented in Sec. II V C] could be obtained from 
more elementary reasoning about uncertainty relations 

In related work, Shabani and Jacobs [13, HH found 
the globally optimal (in time) protocol for D = 3 and 
the locally optimal protocol for all dimensions for reduc- 
tion of a quantity related to impurity. Their calculated 
bounds on the speed-up were 2{D — 1) < Ssj < 2{D—1)'^. 
Naively comparing their bounds to ours suggest that an 
unbiassed protocol is not optimal. However it is not yet 
clear if their upper bound is achievable for D > A. 

It is possible to perform a simple calculation, indepen- 
dent of the one presented in Ref. [53|, for the absolute 
upper bound on the speed-up. First we calculate a lower 
bound on Efdi]: 



E[dL] = -87dt|Tr 



Tr 



pXpX 
pX 



2Tr 



Xp 



Tr 



Xp' 



(64) 



by using the trace inequality Tr [ABAB] < Tr [A^B^] , 
for A, B Hermitan. The breve on X denotes the uni- 
tary transformation X = U^XU, where U is any uni- 
tary. Note that all but the first term, Tr [XpXp], are are 
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0(1 — LP') [or 0(L^)]. Consequently the first term domi- 
nates E[(iL] as it is 0{L). Next we make a two eigenvalue 
approximation to p. We find 

E[dL2] = -8jdtTT\p2Xp2x] (65) 
< -8-/dtTT [x'^UplU^] < 2{D ~ 1)2. (66) 

This implies an upper bound on the speed-up of S = 
2{D — 1)2, as found by Shabani and Jacobs [50] ■ The 
results of a numerical search, like the one performed in 
Fig. [lH indicates that this bound is indeed achievable. 
Furthermore we can confidently say that the speed-up of 
the time-optimal control strategy for impurity reduction 
is bounded above by 5* = 2{D — 1)^. This is justified by 
making a two eigenvalue approximation to p2 and then 
using^ the proof that Jacobs' protocol 3 is optimal [Tl|j 

It is interesting that the rapid measurement proto- 
col (considered in Refs. 0|) and the complementary 
measurement protocols both afford at most a speed-up 
0{D^). They work in very different ways. The rapid 
measurement protocol is essentially classical in nature. It 
uses operations that ensure the state and measurement 
commute at all times. In particular the operations are 
permutations of the state in the measurement basis. The 
phase space picture for the rapid measurement protocol, 
in the long time limit (i.e. under the two eigenvalue ap- 
proximation), would be two Jz eigenstates (Dicke states) 
placed at ±J. As the measurements are in the z ba- 
sis, the mechanism for the rapid measurement speed-up 
is an increase in signal to noise ratio in the measure- 
ment record. In this article we have shown that the com- 
plementary measurement protocol is essentially quantum 
mechanical in nature. The mechanism for the purifica- 
tion and speed-up is measurement backaction (at least for 
the QFT protocol). The two approaches also have differ- 
ent advantages. Rapid measurement has the advantage 
that it enables one to obtain information about the ini- 
tial state. The complementary protocol presented here 
has the advantage that it provides a nearly determinis- 
tic improvement in information gathering of the current 
system state. In future work we will explore if it is pos- 
sible to obtain information about the initial state for a 



continuous complementary measurement. 

Finally we will speak to the practicality of implement- 
ing QFT based protocols. In recent work we showed how 
to analyse imperfections in a purification protocol using 
a feedback master equation [53| so we will not discuss 
such issues further here. Instead we will focus on the 
implementation the QFT. For qbits the QFT is a funda- 
mental quantum logic gate known as the Hadamard gate 
[33} , and is typically easy to implement. 

For D > 2 implementing the QFT is more difficult, as 
it involves multiple logic gates (for n-qbit systems) [5J| or 
multi-level coherent operations (for a qudit). Earlier we 
mentioned that in some atomic spin systems it should 
be possible to construct the desired QFT unitary 41]. 
We think, however, that solid state systems are the most 
likely candidate for which it would be useful to imple- 
ment the ideas presented in this article. This is because 
in solid state systems the measurement strength 7 is typ- 
ically much smaller (by an order of magnitude or more) 
than the maximum control strength |a| (e.g. in Ref. [55| 
1/7 is on the order of microseconds and l/|a| is of order 
nanoseconds). Consequently it is possible to imagine ap- 
plying feedback continuously through the measurement 
process so that the eigenbasis of the state and observable 
are QFT pairs. Alternatively, if the maximum control 
strength is large enough, the feedback could be applied 
impulsively at discrete times, which can gives results sur- 
prisingly close to the continuous version (as noted in 
Fig. Eland Ref. 0). Further the exciting results in the 
field of superconducting qdits (D S [3,5]) where measure- 
ment and full unitary control has been demonstrated for 
a number of systems [5^ [13] lead us to speculate that 
our proposals could be experimentally tested within the 
next decade. 
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